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Abstract 

We study the evolution of the distribution of eigenvalues of a N x N ma- 
trix ensemble subject to a change of variances of its matrix elements. Our 
results indicate that the evolution of the probability density is governed by 
a Fokker-Planck equation similar to the one governing the time-evolution of 
the particle-distribution in Wigner-Dyson gas, with relative variances now 
playing the role of time. This is also similar to the Fokker-Planck equation 
for the distribution of eigenvalues of a N x N matrix subject to a random 
perturbation taken from the standard Gaussian ensembles with perturbation 
strength as the "time" variable. This equivalence alongwith the already known 
correlations of standard Gaussian ensembles can therefore help us to obtain 
the correlations for various physically-significant cases modeled by random 

banded Gaussian ensembles. 
PACS numbers: 05.45+b, 03.65 sq, 05.40+j 
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It is now well-established that the statistics of energy levels of complex quantum systems 
e.g. chaotic systems, disordered systems can be well-modeled by the ensembles of large 
random matrices The nature of underlying quantum dynamics of these systems divides 
the applicable random matrix ensembles (RME) into two major classes. The ensembles of 
the full random matrices (RFME) f|| have been very successful in modelling the statistical 
properties of systems with delocalized quantum dynamics. On the other hand, presence 
of the localized quantum dynamics requires a consideration of the Band random matrices 
ensembles (RBME) which can generally be described as a N x iV matrix with non-zero 
elements effectively within a band of the width b > 1 around the main diagonal. The presence 
of these two major classes can be understood as follows. In the matrix representation of 
an operator, an off-diagonal matrix element describes the overlapping of the eigenfunctions 
in the basis space. For the delocalized quantum dynamics, the eigenfunctions are extended 
and the overlapping between them and therefore off-diagonal matrix elements are of the 
same order. The ensembles of such matrices are referred here as the RFME. In contrast, 
the localization of the eigenfunctions implies a decaying tendency of overlapping between 
them, leading to a diminishing strength of the corresponding off-diagonal matrix elements 
and therefore RBM. Our aim in this paper is to suggest a general method for the statistical 
studies of RBME. 

Various studies of the eigenfunctions of complex systems have revealed the existence of 
various types of localizations, intermediate between fully localized eigenfunctions (occurs 
only when system size L — > oo) and fully extended eigenstates (implying that the local- 
ization length > the system size). For example, the eigenfunctions, in a weakly disordered 
potential exhibiting the Anderson metal-insulator transition, are essentially structureless and 
extended and the statistics here can be well-modeled by FRME. In the insulator regime, a 
complete localization of eigenfunctions leads to a zero correlation among them and a Pois- 
son statistics for energy levels. However the eigenfunctions in the critical region near the 
Anderson transition reveal a special feature "multifractality" in their structure ||. The mul- 
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tifractal nature of a wavefunction is commonly accepted to be related to its space structure 
which also reveals itself in their overlapping with each other; the two fractal wavef unctions, 
however sparse they are, overlap strongly in contrast to two fully localized states. The na- 
ture of the energy level statistics being related to the structure of eigenfunctions as well as 
their overlapping, the latter's special feature should manifest also in the former. Indeed the 
presence of multifractality leads to a new type of universal critical level statistics which is 
described by a set of critical exponents and is different from both Wigner-Dyson (FRME) as 
well as Poisson statistics jH|. Recently there have been suggestions regarding the existence 
of such eigenfunctions in wide variety of physical systems. For example, a Coulomb impurity 
inside an integrable square billiard leads to multifractal eigenfunctions in momentum repre- 
sentation however small is the strength of the potential 0. It also seems quite possible that 
the multifractality is a quantum mechanical manifestation of the special nature of under- 
lying classical dynamics, intermediate between integrability (fully localized eigenfunctions) 
and chaos (fully extended eigenfunctions) [[J . 

Recent research on complex systems has indicated that RBME, with zero mean value of 
all the matrix elements and the variance < \Hij\ 2 >oc a(\i—j\), can serve as good models for 
systems with multifractal eigenfunctions here the function a(r) decays with r either 

as a power law or exponentially (or faster) at r >> 1. However various other RBM type 



structures appear in many other physical contexts, for example, nuclei, atoms [10], solid 



state [|Tl|], quantum chaos |12[]. One such example is that of WRBM (W stands for Wigner) 
|| where the mean value of diagonal elements increases linearly along the main diagonal 
(< H nn >= an). These matrices have been shown to be a good model, for example, for 
the tight-binding Hamiltonian of a quantum particle in a 1-D disordered system subject to 



a constant electric field as well as for heavy atoms and nuclei [[0]. Another example is of 
the ensemble of banded matrices with diagonal elements fluctuating much stronger than the 
off-diagonal ones (< \Ha\ 2 > / < \Hij\ 2 >oc b » 1 ), also known as RBM with preferential 
basis. A very interesting problem of two interacting particles propagating in a quenched 
random potential can effectively be mapped on to this class Jl4| . 



The nature of a localized dynamics leaves its imprints on the distribution properties of 
the matrix elements of the generator of motion e.g. Hamiltonian. It is therefore desirable to 
study RBME with various distributions and recently a lot of effort has been applied in this 
direction. However the techniques suggested so far [|^[TT]|L3| depend strongly on the type of 
distribution chosen and mostly give approximate results. Our aim in this paper is to suggest 
a method which leads to results in a generic form vaild for most of the RBM models. This we 
achieve by mapping the problem to the study of a particular class of transition ensembles, 
the one arising during the Poisson — > standard Gaussian transitions under a perturbation 
taken from the standard Gaussian ensembles (SGE). The pre-existing information regarding 
the correlations in the latter case therefore can help us to obtain the same for RBME. 

We proceed as follows. Our interest is in the evolution of the eigenvalue distribution 
of H, taken from an ensemble of hermitian matrices, due to variation of the variances 
of its matrix elements. We choose the distribution p(H) of matrix if to be a Gaussian, 
p(H,y) = Cexp(—J2k<i a kiHl l ) with C as the normalization constant and y as the matrix 
of relative variances y k i = °^g£l . Let P({/ij}, y\H ) be the probability of finding eigenvalues 
Aj of H between pi and pi + dpi at a given y (with H as an initial condition), 

N 

P({^}, y\H ) = I J] Km - Xi)p(H, y)dH (1) 

As the a-dependence of P in eq.(l) enters only through p{H) ( — 
[(2a fci ) _1 — Hlj\p{H)), this equality followed by a repeated use of the partial integration 
alongwith eigenvalue equation H = T AO, with A as the eigenvalue matrix and O an 
orthogonal matrix, leads to following, 

>dP({in},y) 

*H ' 1 kl 

where 



2 = —hi (2) 
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Ik, -EtT- / dH P( H )l^r- [WtM-^OnhOm] (3) 
n op n J on k i \ i j 

Let us first study the case with same variance for all the diagonal matrix elements such 

that guoiki = ot for k = I (g k i = 1 + 5 k i) while keeping the variances (=a^ 1 ) of the off- 



diagonals arbitrary. By expressing J2k<i - - = z £fc<J Ek<i and using 

eq.(2), we obtain the following relation 

dP 1 1 
2 J2( 1 ~Vki)Vki^ — = Y. J M Hhi (4) 

k<i d Vki k <i 9kiotu a k ^ 

The left hand side of above equation, summing only over all distinct an, can be rewritten 
as |^ with Y given by the condition that = 2J2k<iVki(^- — Uki)^^- The first term on 
the right hand side (eq.(4)) can further be simplified by first using a partial integration and 
subsequently the eigenvalue equation for H: J2k<i z~Jki = J2n jfp~ i^nP)- The second term 
can similarly be reduced by differentiating the terms inside the brackets in (3) giving us two 
integrals. The use of orthogonality relation of matrix O in the first integral so obtained and 
the equality J2k<i d °dH° nl = J2m A \ x in the second |l5j leads to a F-P equation 
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a v d/i r , 



9 + £ ; 



(5) 



dfi n m ^ n fi n fl r , 

where (3=1. By using the unitarity of eigenvectors and following the same steps, it can 
be proved for complex Hermitian case too (now (3 = 2). It should be noted here, for the later 
use, that the required form of eq.(5) is obtained by applying an appropriate partitioning of 
sums appearing in eq.(4) which leads to rewriting of the terms containing all the unequal 
coefficients a^i as the derivative with respect to parameter Y . The coefficient appearing 
with drift and diffusion terms in eq.(5) is the one common to many of the matrix elements. 

The definition of Y depends on the relative value of variances of the off-diagonals as 
discussed below for some cases: 

Case I. When all the off-diagonals have same variances such that a^i = a' (for k ^ I) 

In this case, as y — ^, therefore Y = 2 log-^-q-. 

Case II. When variances of the off'- diagonals change with respect to distance from the 
diagonal such that an = ct T (for r = \k — l\ > b > 0) and otki — ct for r < b. 

Now all the off-diagonals (with y r ^ 1) contribute separately to Y, Y = 2 Y^^ =b+1 log i| 
where y r = — . 



For all y r > 1, this represents a RBM with strongly fluctuating diagonal elements (known 
as RBM with preferential basis). The case with all y r < 1 describes the standard RBM, 
used to model, for example, the spectral statistics of the tight binding Hamiltonian of a 
quantum particle in a 1-D system with long range random hoppings and also for quasi 
1-D disordered wires JTE|] . 

Case III. When the off- diagonals form various groups with different varainces (with same 
variance for each matrix element in one such group). 

Now contribution to Y comes from each such group. For example, with M groups of 
variances given by a" 1 (g — 1 — > M, a g ^ 1), Y — 2 Z^li log ^^ where y g = This will 
be helpful to model the more general cases where the variance < \Vij\ 2 > is dependent on 
both indices i,j instead of their difference \i — j\. 

The steady state of eq.(5), P({/ij}, oo) = P^ = fli<j — ¥ ^ fcM *, is achieved 

for Y — > oo which corresponds to y — > 1 in case (I), y r — * 1 in case (II) and yu — > 1 
for the case (III). This indicates that, in the steady state limit, system tends to belong to 
the standard Gaussian ensembles. Eq.(5) (later referred as variance- variation or VV case) 
is formally the same as the F-P equation governing the Brownian motion of particles in 
Wigner-Dyon gas [[J with transition parameter being the relative variance in the former and 
time in the latter. This is also similar to the F-P equation for the eigenvalue distribution 
of a hermitain matrix H — H Q + rV undergoing random perturbation V of strength r 
and taken from a SGE with arbitrary initial condition H (later referred as perturbation 
variation or PV case) [JD]]; here r acts as the transition parameter. P({/ij}, Y) can therefore 
be obtained by the same procedure as used in PV case which is briefly given as follows. The 



transformation \& = -P/\APoo reduces eq.(5) to a "Schrodinger equation" form, |p = 



where the 'Hamiltonian' H turns out to be a Calogeo-Moser (CM) Hamiltonian |T9|,|20 



H = Y,i — \ J2i<j (f^rja — x^*/^' an d nas well-defined eigenstates and eigenvalues 
[PHI . The "state" \I/ or P({fii}, Y\H ) can therefore be expressed as a sum over its eigenvalues 
and eigenfunctions which on integration over all the initial conditions Ho leads to the joint 
probability distribution P({/ij}, Y) and thereby correlations. Unfortunately, due to technical 



problems, the latter could be evaluated only for the transitions with final steady state (limit 
A — > oo) as GUE p^j . However a set of hierarchic relations among the correlators for all 
transitions in PV case, and therefore for VV case, can be obtained by a direct integration 
of the F-P equation for P (eq.(5)) [P!5| , p0[| . For example, in large iV-limit, the evolution 
of level density p(/i, Y) is governed by the Dyson-Pastur equation [|T§|,pU| which results 



in a semi-circular form for p (thus agreeing well with the p obtained in [|18j by super- 
symmetry technique). The relation p(p,Y) = N~ 1 Y^ n p(fjL,n,Y) can further be used to 
obtain the equation for the local density of states (LDOS) p(p,n,Y), a more informative 
and experimentally accessible quantity (referred below as p n ), 

§^-4(? p / d ^K < e) 

This nonlinear equation gives uniquely p(E, n, Y) starting from an initial p(E, n, 0). 

For localization studies, it is appropriate to choose the initial ensemble Hq as that of 
diagonal matrices (P(H ) oc e~2Yui H u with Va = p 0i and all y — > oo) which corresponds 
to Poisson distribution for eigenvalues and Y = 0. With equilibrium distribution (Y — > oo) 
given by SGE, this case thus represents a Poisson — > SG transition with Y as a transition 
parameter and the intermediate ensembles representing various RBMEs depending on the 
type of y r 's. Thus the correlations here will be similar to those in the Poisson— >SG transition 
in SGE (PV case with H taken from a Poisson ensemble). We already know that the 
transition for PV case is abrupt for large dimensions and finite r and a rescaling of r by 
mean spacing D 2 ((x N) is required to make it smooth, the new transition parameter being 
A = -jp. A similar rescaling should also be applied to Y in VV case for the same reason; 
here too D 2 oc N |18|]. (This A^-dependence of D also follows from analogy of the evolution- 
equations for level-densities in the two cases ). 

Fortunately the two-point correlation R2{r; A) for Poisson — > GUE transition in PV case 
has already been obtained [^(J and, as discussed above, is also valid for the VV case (now 
A = Y/D 2 ) . 
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^ 2 )(l + 2^ v /i)' 



i? 2 (r; A) — Ri(t\ oo) = — dx dz cos(2nrx) exp — 87r 2 Aa;(l + x + 2z^fx) . 

n Jo J-i L J \ l+x + 2z\/x 

where R2(r,oo) 

= i _ E^l (the GUE limit). As can easily be checked, above equation 
has the correct limiting behaviour, that is, R2 = for A — ► (the Poisson case Y — 0) 
and i?2 = -^2(7"; 00) for A — > 00. As obvious from eq.(7), R2 for intermediate ensembles 
will depend on definition of Y and therefore on the nature of localization which results in 
various types of level-statistics. For example, let us calculate the correlation for one such 
case, namely, Hij = Gija(\i — j\) with G a typical member of SGE and a(r) = 1 and (b/r) a 
for r < b and > b (b » 1) respectively. For this case, y r = (|) CT and Y ps 2£]^U + i (r) 
which gives A oc N 2 ~ 2a . For o > 1 with — > 00, therefore, the eigenvalue statistics 
approches Poisson limit, A being very small. Similarly for a < 1, A is sufficiently large and 
the eigenvalue statistics approches SG limit. For a — 1, the independence of A from A^ leads 
to an eigenvalue statistics, very different from that of SGE or Poisson and therefore agrees 
well with the results obtained in by using non-linear a-model technique. Further R 2 , 
(eq.(7)) is also in good accordance with the one given in [0] (which can be seen by a direct 
substitution of eq.(52) of in eq.(17) of 0). 

It is possible to have physical situations when the matrix elements form various groups 
such that those in one group have the same value for akign = « r (r = 1 — ► M with M as total 
number of groups). This case can similarly be treated by using 2J2rLi — a r^~ = 
-^-J2 r R — J2r with I r = J2{k,i}e r Iki and Iki still given by eq.(3). Here a s refers to the 
a-coefficient of any one (chosen arbitrarily) of the M groups. This partitioning again leads 
to eq.(5) (with a — >• a s on the right hand side) and therefore similar level-correlations in 
terms of the transition parameter Y = 2 logflr \ y V -i\ with y r — ^r- 

Finally let us consider the case of WRBM when mean value of the diagonal elements 
increases along the main diagonal: < H nn >= 7/(n). For simplification, we still take p{H) 
to be a Gaussian: p{H) = Ce~^- Jk < l ° lkll ' Hkl ~ 1 ^ k ^ &k ^ with all distinct a^. Proceeding as 
before, one again obtains eq.(2) with Ri given by eq.(3) but now, instead of applying the 
partitioning of the sum, we just evaluate 2Y,k<i ~^~ l a ti9ki- This leads to a F-P equation 



without the linear drift term, 



dY 



E 



n 



d\i n <9/i, 



d d 



n 



P 



(8) 



m- 



where Y = J^kKii^kigki)' 1 and the steady state P^ = Yli<j 1/^ ~ Vjf is obtained in limit 
Y —>■ oo or aki — > 0. Similarly for the degenerate case, when all the au can be divided 
into M groups with distinct values a r = akign, r = 1 — > M, 2Y,rLi §^~ a r can be used 
to obtain eq.(8) where now Y = Y^ r (2a~ ). An application of the same transformation 
\& = P/y/Poo reduces eq.(8) in the " Schroedinger equation" form with Hamiltonian H, 



confining potential and, as shown in |[20|| , has essentially the same solution for \1/ (without 
the Gaussian factor) and therefore P(fi, Y\H$). However, for the evaluation of Y), one 
needs to integrate over all the initial conditions through which the non-zero mean values 
of the diagonals enter in the calculation and may influence the statistics for large 7-values. 
Nontheless, in large iV-limit, the hierarchic relations among the correlators are of the same 
type as in other RBM cases [0,^] with LDOS satisfying the eq.(6) |P0[j . Here the choice of 
H from a Poisson ensemble (p(H ) oc e'^Y.j^n-'yfU)) ^ gj ves the initial value of Y = a -1 . 
With Y — > 00 corresponding to a SG type ensemble, therefore, WRBM also appear as the 
intermediate ensembles in the Poisson — > SG type transition. 

In this paper, we have analytically studied the response of energy levels of complex 
quantum systems to the changing distribution of matrix elements of their hamiltonians. 
Our results indicate that the n th order correlations at a given variance are the same as the 
corresponding eigenvalue-correlations of hamiltonian H = Hq + tV at a given r value with 
H taken from a Poisson ensemble and V taken from a SG ensembles. This analogy also ex- 
tends to the 2 nd order parametric correlators however the Markovian nature of F-P equation 
restricts from making similar conclusions about higher orders. The intermediate ensembles 
arising in the VV transition correspond to various types of RMEs, and as our mapping 
suggests, the results for most of them can be obtained by studying just one transition in full 
detail, namely, an appropriate initial ensemble — > SGE caused due to a SG perturbation; the 



d 2 



— \J2i<j ( _ which is same as CM Hamiltonian except for the absence of 
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knowledge of their correlations can help us in statistical studies of many important physical 
properties of complex systems. Further, as discussed in || , some of the RBMEs can also be 
connected to various other types of ensembles; a knowledge of the properties of former will 
therefore help in statistical analysis of the latter. It should also be possible to apply this 
technique to the ensembles of Unitary matrices (periodic RBMs) as well. 

Our study still leaves many questions unanswered. At present, it is not clear whether our 
method can be extended to non-Gaussian distributions as well, at least in some limit. We are 
also unable to say whether a similar mapping can also be done for the eigenvector statistics. 
In this connection, however, it should be mentioned that analogy of equations governing the 
eigenvalue distribution leads to a similar form of equations for the distribution of matrix 
elements @ and, therefore, it should somehow manifest itself in eigenvector distributions 
too. 

I am greatful to B.S.Shastry, B.L.Altshuler and V.Kravtsov for useful suggestions during 
the course of this study. 
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